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Abstract 

Modeling  and  testing  of  inflatable  membrane  structures  for  rapidly  deployable  port  infrastructures  is 
of  interest  to  the  US  Navy’s  new  capabilities  development.  For  decades,  the  Navy  has  used  small 
displacement  landing  craft  to  be  the  conduit  between  heavy-tonnage  ships  offshore  and  the  landing 
zone  ashore.  This  has  been  a  means  of  transporting  troops  and  their  combat  and  logistics  support.  A 
problem  with  this  procedure  is  that  displacement  landing  craft  can  only  get  close  enough  to  the 
beach  to  disembark  vehicles  and  personnel  to  wade  the  final  distance  ashore.  Although  the  water  is 
shallow,  the  vehicles  need  to  be  designed  to  operate  in  a  corrosive  seawater  environment.  An 
alternative  method  the  Navy  considers  is  floating  causeways. 

The  major  objective  of  this  study  is  to  assess  the  feasibility  of  using  an  advanced  numerical  code  to 
simulate  the  dynamics  of  a  mobile  support  platform  for  the  deployment,  retrieval,  and  operation  of 
the  MOSES  inflatable  causeway  under  wave  action.  This  is  achieved  by  developing  a  numerical 
model  of  the  wave  basin  with  wave  generation  over  a  rigid  basin  bottom  to  study  the  dynamic 
behavior  of  MOSES  under  laboratory  conditions.  The  purpose  of  this  study  is  to  compare 
qualitatively  the  numerical  predictions  with  experimental  results  obtained  in  the  experimental  study 
to  assess  the  feasibility  of  such  modeling. 
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Executive  Summary 


MOSES  is  a  conceptual  inflatable  causeway  system  designed  to  provide  access  for  shallow  draft 
vessels  near  the  shore  to  the  beach  and  it  is  designed  to  be  lightweight  and  easy  to  deploy.  The  large 
fabric  membrane  structure  is  filled  with  seawater  and  pressurized  to  become  relatively  rigid.  This 
causeway  rests  on  the  seafloor  forming  a  roadway  about  a  meter  above  the  sea  surface. 

The  major  objective  of  this  study  is  to  assess  the  feasibility  of  using  an  advanced  numerical  code  to 
simulate  the  dynamic  motion  of  MOSES  under  wave  action.  This  is  achieved  by  developing  a 
numerical  model  of  the  wave  basin  with  wave  generation  over  a  rigid  basin  bottom  to  study  the 
dynamic  behavior  of  MOSES  under  laboratory  conditions.  The  numerical  predictions  were 
compared  qualitatively  with  experimental  results  obtained  in  an  accompanying  experimental  study. 
The  simulations  thus  far  utilize  an  Arbitrary  Lagrangian  and  Eulerian  (ALE)  technique  to  capture 
the  multi-physics  phenomenon. 

MOSES  was  treated  as  a  fully  rigid  object  and  then  a  single  flexible  case  was  tested  to  observe  its 
motion  under  the  action  of  the  waves.  Results  show  that  the  coupled  FSI  numerical  code  can 
successfully  capture  the  motion  of  MOSES  under  the  severe  action  of  waves.  All  the  numerical  test 
cases  described  above  show  good  matching  of  the  physics  of  the  wave  motion  and  the  motion  of  the 
membrane  structure  on  the  movable  beach  and  the  qualitative  behavior  compared  well  with  those 
shown  in  the  experimental  test  cases.  Both  the  numerical  test  cases  of  treating  the  water-filled 
membrane  (MOSES)  as  rigid  and  fully  flexible  show  a  realistic  behavior  of  the  membrane  structure 
motion  under  the  action  of  waves. 

A  detailed  descriptions  of  the  numerical  simulation  procedure  specific  to  the  modeling  the  physical 
wave  basin  and  water-filled  membrane  (MOSES)  using  a  multi-physics  multi-scale  finite-element 
code  is  also  delineated.  Specifically,  a  nine-step  simulation  procedure  pertaining  to  such  class  of 
fluid-structure  interaction  problems  are  shown. 
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INTRODUCTION 


Modeling  and  testing  of  inflatable  membrane  structures  for  rapidly  deployable  port 
infrastructures  is  of  interest  to  the  US  Navy’s  new  capabilities  development.  For  decades,  the 
Navy  has  used  small  displacement  landing  craft  to  be  the  conduit  between  heavy-tonnage  ships 
offshore  and  the  landing  zone  ashore.  This  has  been  a  means  of  transporting  troops  and  their 
combat  and  logistics  support.  A  problem  with  this  procedure  is  that  displacement  landing  craft 
can  only  get  close  enough  to  the  beach  to  disembark  vehicles  and  personnel  to  wade  the  final 
distance  ashore.  Although  the  water  is  shallow,  the  vehicles  need  to  be  designed  to  operate  in  a 
corrosive  seawater  environment.  An  alternative  method  the  Navy  considers  is  floating 
causeways. 


Physical  Modeling 


MOSES  is  one  such  conceptual  inflatable  causeway  system  designed  to  provide  access  for 
shallow  draft  vessels  near  the  shore  to  the  beach.  The  system  is  designed  to  be  lightweight  and 
easy  to  deploy  (Bloxom  et  al  2010).  MOSES  is  a  large  fabric  membrane  structure  filled  with 
seawater  and  pressurized  to  become  relatively  rigid.  This  causeway  rests  on  the  seafloor  forming 
a  roadway  about  a  meter  above  the  sea  surface.  The  goal  of  the  MOSES  project  is  to  design  a 
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causeway  that  can  reach  large  ships  several  hundred  meters  from  the  shore,  keep  vehicles  dry, 

and  be  operational  in  at  least  Sea  State  4  (SS4)  conditions.  On-going  research  of  inflatable 
structures  shows  that  they  can  be  easily  deployed,  lightweight  and  compact  making  them 
appealing. 


Using  the  seafloor  for  support  enables  the  causeway  to  remain  stable  in  high  sea  states  remaining 
relatively  unaffected  by  large  swells.  The  weight  of  the  water  in  the  membrane  structure  above 
sea  level  creates  a  large  downward  force,  which  holds  the  system  in  place.  To  minimize  rolling, 
MOSES  contains  internal  fabric  ribs  that  run  its  entire  length.  These  ribs  create  a  structure  with  a 
flat  bottom  and  top.  The  selected  shape  provides  for  a  roadway  on  top  and  the  flat  bottom 
ensures  stability  against  rolling.  An  experiment  on  a  scaled  model  of  MOSES  was  conducted  in 
the  summer  of  2010  (Bloxom  et  al  2010). 


Numerical  Simulation 

Large-scale  Naval  fluid  structure  contact-impact  problems  such  as  the  conceptual  inflatable 
causeway  project  form  an  interesting  and  challenging  interdisciplinary  area  of  Fluid-Structure 
interaction  (FSI).  Methods  that  use  results  from  physical  wave  basin  experiments  to  assess  the 
wave  forces  on  such  structures  are  usually  employed  to  investigate  such  class  of  problems.  These 
experiments,  though  closer  to  the  physical  phenomena,  often  tend  to  be  time-consuming  and 
expensive  and  it  is  not  economically  feasible  to  conduct  parametric  studies.  Alternatively, 
numerical  models  when  developed  with  similar  capabilities  will  complement  the  experiments 
well  because  of  the  lower  costs  and  the  ability  to  study  phenomena  that  are  not  feasible  in  the 
laboratory.  LS-DYNA  is  one  such  numerical  model  that  is  used  most  often  to  analyze  this  class 
of  contact-impact  problems. 

LS-DYNA  (developed  by  LSTC)  is  a  multi-functional  code  with  explicit  and  implicit  finite- 
elements  to  simulate  and  analyze  highly  nonlinear  physical  phenomenon  which  usually  involves 
large  deformations  within  short  time  duration.  It  is  an  advanced  nonlinear  finite-element  based 
commercially  available  numerical  software  which  has  arguably  the  best  solid  mechanics,  contact 
and  impact  models  among  commercially  available  software,  and  has  growing  fluid  mechanic 
features  (Hallquist,  2006).  The  main  solution  methodology  is  based  on  explicit  time  integration. 
An  implicit  solver  is  currently  available  with  somewhat  limited  capabilities  including  structural 
analysis  and  heat  transfer.  It  also  has  an  excellent  computational  framework  combined  with  pre- 
and  post-processing  capabilities  to  address  various  research  and  application  needs.  A  significant 
number  of  fluid-structure  interaction  models  in  the  industry  have  been  successfully  employed 
using  finite-element  (FE)  codes  with  an  arbitrary  Lagrangian-Eulerian  (ALE)  formulation  and  a 
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contact-impact  algorithm  because  this  code  enables  the  modeling  of  the  instantaneous  coupling 

between  the  fluid  and  the  structure  (Hallquist,  2006). 


Contact-impact  algorithms  have  always  been  an  important  stand  out  capability  in  the  LS-DYNA 
family  of  codes.  Contact  and  impact  may  occur  along  surfaces  of  a  single  body  undergoing  large 
deformations,  between  two  or  more  deformable  bodies,  or  between  a  deformable  body  and  a 
rigid  barrier.  This  type  of  interface/contact  is  mainly  used  when  there  is  a  need  for  the  ALE 
transformations  for  the  impact  scenario  (specifically  meant  for  FSI  problems).  By  a 
specialization  of  this  algorithm,  such  interfaces  can  be  rigidly  tied  to  admit  variable  zoning 
without  the  need  of  mesh  transition  regions  (Hallquist,  2006).  The  purpose  of  the  present  study  is 
to  provide  a  preliminary  assessment  of  an  advanced  numerical  code  in  simulating  complex  fluid- 
structure  interaction  behaviors  of  the  model-scale  system. 


Objectives 


The  major  objective  of  this  study  is  to  assess  the  feasibility  of  using  an  advanced  numerical  code 
to  simulate  the  dynamics  of  a  mobile  support  platform  for  the  deployment,  retrieval,  and 
operation  of  the  MOSES  inflatable  causeway  under  wave  action.  This  is  achieved  by  developing 
a  numerical  model  of  the  wave  basin  with  wave  generation  over  a  rigid  basin  bottom  to  study  the 
dynamic  behavior  of  MOSES  under  laboratory  conditions.  The  purpose  of  this  study  is  to 
compare  qualitatively  the  numerical  predictions  with  experimental  results  obtained  in  the 
experimental  study  (Bloxom  et  al  2010)  to  assess  the  feasibility  of  such  modeling. 


SUMMARY  OF  PHYSICAL  MODEL  AND  EXPERIMENT 


MOSES  Physical  Model 


MOSES  is  designed  to  be  pressurized  to  make  the  system  sufficiently  rigid  for  vehicles  to  drive 
across  (roll-on  roll-off).  The  required  pressure  is  induced  through  the  use  of  reservoirs  that  are 
open  to  the  atmosphere.  The  reservoirs  serve  two  main  purposes.  Firstly,  they  raise  the  level  of 
the  water  above  the  roadway,  which  creates  a  sufficient  pressure  head  to  make  the  MOSES  rigid. 
Secondly,  the  reservoirs  act  as  walls,  deflecting  the  large  swells  expected  in  SS4,  a  design 
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requirement.  A  secondary  benefit  of  the  reservoirs  is  that  they  provide  some  allowance  for 

punctures  or  tears  in  the  fabric.  The  large  reservoirs  will  enable  the  causeway  to  remain 
sufficiently  rigid  for  several  minutes  after  a  significant  puncture  to  allow  time  for  safe  evacuation 
of  vehicles  and  personnel. 


Since  the  system  needs  to  be  easily  deployed,  a  simple  mechanism  has  been  designed  to  unroll 
the  deflated  causeway  from  the  transport  ship  to  the  shore.  Using  fabric  as  the  main  structural 
component  allows  MOSES  to  be  easily  compacted  and  rolled  around  a  giant  spool.  The  recovery 
process  is  simply  the  reverse  of  the  deployment  process. 


The  viability  of  the  MOSES  causeway  concept  was  tested  using  a  scale  model.  The  MOSES 
system  is  novel  because  of  its  fabric  material  construction,  which  can  be  compacted,  stowed  and 
easily  deployed.  Compared  to  the  current  causeways,  MOSES  has  advantages  in  terms  of 
deployment/recovery  times,  stowage  volumes,  and  overall  weight.  These  attributes  make  the 
inflatable  causeway  appealing  for  future  development  of  this  ship  to  shore  transition  technology. 
(Dickens  and  Rosen,  2007) 


Experimental  Facility 


As  described  in  Bloxom  et  al  (2010)  an  experimental  investigation  of  the  fluid-structure 
interaction  of  water  filled  inflatable  membrane  structure  in  the  near  shore  environment  was 
performed  by  undergraduate  students  from  a  number  of  universities  across  United  States 
supported  by  a  naval  summer  research  program  (Figure  1).  The  structure  of  interest  was  a  10'  x  2' 
x  0.75'  tubular  bag  developed  at  the  Center  for  Innovation  in  Ship  Design  (C1SD)  at  the  Naval 
Surface  Warfare  Center  (NSWC),  Carderock  Division  as  a  proof  of  concept  for  the  design  of  a 
rapidly  deployable  inflatable  structure  causeway  to  be  used  either  as  a  ship  to  shore  connector  or 
a  breakwater. 

The  experiments  were  performed  over  a  range  of  test  conditions  including  three  incident  wave 
angles,  three  water  depths,  and  a  number  of  wave  heights  corresponding  to  various  sea  states. 
Results  confirmed  the  previous  conclusions  that  the  bag  is  stable  and  well  grounded  for  most 
operational  sea  conditions.  Large  amplitude  and  low  frequency  waves  can  induce  significant 
motions  of  the  structure,  but  the  static  and  dynamic  frictional  coefficients  between  the  structure 
and  the  surface  in  contact  play  a  critical  role  in  these  motions.  For  conditions  where  the  structure 
was  at  an  angle  of  45°  to  the  incident  waves,  highly  nonlinear  wave  conditions  are  produced 
which  created  wave  over-topping  and  oscillatory  motions  of  the  structure. 

The  experimental  test  configuration  was  classified  in  two  different  categories.  One  that  involves 
a  beach  set  up  and  the  other  in  the  absence  of  a  beach.  Figure  2  shows  the  test  configuration  in 
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the  presence  of  a  beach  for  different  orientations  of  the  beach.  Similarly,  Figure  3  shows  the  test 

configuration  in  the  absence  of  the  beach  with  various  degrees  of  orientation.  The  dimensions  of 
the  USNA  wave  basin  and  the  details  of  the  wavemaker  are  summarized  as  follows: 


•  Practical  water-depth  range 

-  6”  -  14” 

•  Movable  beach 

-  Height -9.75” 

-  Width -96” 

-  Length  -  82” 

-  Slope- 1:8.4 

•  Piston  type  wavemaker 

-  6”  stroke  max. 

-  2  Hz  max. 

-  Regular  waves 

-  Simple  irregular  wave  spectra 
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23' 


Figure  1  Wave  basin  dimensions 
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The  test  configuration  with  the  removable  beach  is  shown  below. 


0  Degrees 


45  Degree 


Figure  2  (a) 


Figure  2  (b) 


90  Degrees 
Figure  2  (c) 

Figures  2(a),  2(b)  and  2(c)  Orientation  of  movable  beach  and  membrane  structure 
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The  test  configuration  without  the  removable  beach  is  shown  below. 


0  Degrees 
Figure  3(a) 


45  Degree 
Figure  3(b) 


90  Degree 
Figure  3(c) 

Figures  3(a),  3(b)  and  3(c)  Orientation  of  membrane  structure  without  beach 
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Experimental  Setup 


The  experimental  and  basin  setup  can  be  summarized  in  two  figures  as  follows.  Figure  4  shows 
the  constant  pressure  reservoir  to  maintain  constant  pressure  head  in  the  membrane  structure. 
Water  was  supplied  steadily  from  one  hose  and  overflow  was  drained  into  another  one  connected 
to  the  bucket.  Figure  5  shows  the  back  view  of  the  wavemaker  used  for  wave  generation  in  the 
experimental  wave  basin  (see  Bloxom  et  al  2010  for  a  detailed  description). 


Figure  4  Constant  pressure  reservoir  to  maintain  constant  pressure  head 


Figure  5  Back  view  of  the  wavemaker  used  for  wave  generation  in  the  wave  basin 
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NUMERICAL  MODELING  AND  SIMULATION  OF  INFLATABLE  CAUSEWAY 


The  complexity  of  the  experiments  demands  an  advanced  robust  predictive  numerical  code  to 
model  the  system  components  accurately.  Numerical  modeling  of  the  MOSES  inflatable 
causeway  is  carried  out  using  a  finite  element  code.  It  is  a  state-of-the-art  nonlinear  multi¬ 
physics  computational  mechanics  code  with  advanced  contact  and  impact  modeling  capabilities. 
The  multi-physics  feature  allows  the  modeling  of  coupled  fluid-structure  interaction  between  the 
waves  and  the  membrane  structure  and  the  contact  and  impact  capability  is  essential  for 
modeling  the  wave  impact  effects.  The  research  group  at  Oregon  State  University  has  a  long¬ 
term  research  collaborative  relationship  with  Livermore  Software  Technology  Corporation 
(LSTC),  the  software  company  that  owns  LS-DYNA. 


Fluid  Models 

Throughout  the  past  decade  considerable  progress  has  been  made  in  implementing  various  fluid 
models  in  the  numerical  code.  Importantly,  the  fluid  flow  solver  capability  is  adopted  to  simulate 
the  viscous  flow  in  this  project.  To  date  these  fluid  models  are  used  for  different  applications, 
varying  from  regular  elastic  and  acoustic  fluid  elements  to  the  complex  null  material  model 
which  demands  an  equation  of  state  to  be  prescribed  (Hallquist,  2006).  A  detailed  description  of 
the  fluid  constitutive  models,  along  with  their  governing  equations  is  provided  in  the  theory 
manual  (Hallquist,  2006).  The  material  model  parameters  that  can  be  calibrated  to  these  data  sets 
have  only  been  used  in  the  existing  models.  Only  those  fluid  models  whose  material 
characterization  is  readily  available  and  widely  tested  for  these  sects  of  problems  have  been  used 
for  the  present  FSI  problem.  The  current  fluid  models  are  based  on  Navier-Stokes  equation,  with 
an  arbitrary  Lagrangian-Eulerian  formulation.  An  equation  of  state  is  essential  to  capture  the 
fluid  behavior. 


NWT  formulations  and  ALE  description 

The  water  domain  in  the  Numerical  Wave  Tank  (NWT)  is  governed  by  the  compressible  Navier- 
Stokes  equations  and  solved  numerically  using  an  ALE  formulation  over  the  entire  domain. 
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(i) 


where  p  is  the  fluid  density,  u  is  the  water  particle  velocity,  p  is  the  bulk  viscosity  and  g  is  the 
body  force  due  to  gravity.  Figure  16  shows  the  computation  domain  of  the  NWT  with  the  FE- 
ALE  mesh  (Appendix). 

Simulations  involving  the  NWT  are  often  demanding  because  of  the  large  distortions  of  the  fluid 
elements  and  there  is  always  a  need  to  track  the  water  (unsteady)  free  surface  precisely.  The 
ALE  formulation  takes  advantage  of  two  classical  descriptions  of  motion:  Lagrangian  and 
Eulerian.  The  Lagrangian  description  is  mainly  employed  in  structural  mechanics  because  of  it 
excellent  capability  in  tracking  free  surfaces  and  interfaces  between  different  materials.  On  the 
other  hand,  the  Eulerian  description  is  widely  used  in  computational  fluid  dynamics  because  it 
can  efficiently  handle  large  distortions  due  to  fluid  particles  movements.  The  ALE  formulation 
can  be  thought  of  as  an  algorithm  consisted  of  a  Lagrangian  time  step  followed  by  a  remap  step. 
The  remap  step  stops  the  calculation  when  the  mesh  is  distorted.  It  then  smoothes  the  mesh,  and 
remaps  the  solution  from  the  distorted  mesh  to  the  smoothed  mesh  (Belytschko  et  al  1982). 


In  theory,  the  Eulerian  description  is  a  subset  in  the  ALE  formulation.  Early  ALE  formulations 
were  not  able  to  smooth  the  mesh  successfully  while  most  of  the  difficulties  have  been 
conquered  by  associating  the  mesh  with  the  nodes  on  the  material  boundaries.  The  ALE 
differential  form  of  the  conservation  equations  for  mass,  momentum,  and  energy  are  readily 
obtained  from  the  corresponding  Eulerian  forms  (Belytschko  et  al  1982): 


Mass: 

Momentum: 

Energy: 


dp  _  dp 
dt  dt 


dv  ,dv 

p —  =  p( — 

dt  dt 


dE  tdE 

p —  =  p( — 

dt  dt 


+  c  ■  S/p  =  -pV  ■  v 

+  (c  •  V)v)  =  V  •  <j  +  pb 
+  c  ■  WE)  =  W -((T-v)  +  v-  pb 


(2) 

(3) 

(4) 
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Where  P  is  the  mass  density,  v  is  the  material  velocity  vector,  <y  denotes  the  Cauchy  stress 

tensor,  b  is  the  specific  body  force  vector,  E  is  the  specific  total  energy,  and  c  is  the  convective 


velocity  c  ~v  v ,  where  v  is  the  grid  velocity. 


The  numerical  wave  tank  ideally  has  a  non-reflective  boundary  condition  at  the  far  end  of  the 
solution  domain  in  order  to  absorb  the  wave  energy  and  prevent  reflection  but  for  the  present 
purpose  a  reflecting-type  boundary  is  used  on  the  far  end  of  the  NWT  to  obtain  the  motion  of  the 
mine  under  a  reflective  wave.  The  other  boundaries  of  the  solution  domain  are  along  the  side 
walls  of  the  tank  (impermeable  no-slip  boundary). 


Structural  Models 


The  code  has  a  wide  range  of  material  and  equation  of  state  models,  each  with  a  unique  number 
of  history  variables.  There  are  approximately  100  material  models  which  can  be  implemented 
according  to  the  complexity  of  the  problem  (Hallquist,  2006).  Spatial  discretization  is  achieved 
by  the  use  of  four  node  tetrahedron  and  eight  node  solid  elements,  two  node  beam  elements, 
three  and  four  node  shell  elements,  eight  node  solid  shell  elements,  truss  elements,  membrane 
elements,  discrete  elements,  and  rigid  bodies.  It  also  has  a  wide  variety  of  element  formulations 
for  each  element  type.  Specialized  capabilities  for  airbags,  sensors,  and  seatbelts  have  tailored 
the  code  for  applications  in  the  automotive  industry.  Currently  the  code  contains  approximately 
one-hundred  constitutive  models  and  ten  equations-of-state  to  cover  a  wide  range  of  material 
behavior.  A  theoretical  manual  has  been  written  to  provide  users  and  potential  users  with  insight 
into  the  mathematical  and  physical  basis  of  the  code  (Hallquist,  2006). 


Contact  and  Impact  Models 


To  enable  flexibility  for  the  user  in  modeling  contact,  the  advanced  numerical  code  presents  a 
number  of  contact  types  and  a  number  of  parameters  that  control  various  aspects  of  the  contact 
and  impact  treatment.  Contact  treatment  forms  an  integral  part  of  many  large-deformation 
problems.  Accurate  modeling  of  contact  interfaces  between  bodies  is  crucial  to  the  prediction 
capability  of  the  finite  element  simulations.  It  also  offers  a  large  number  of  contact  types.  Some 
types  are  for  specific  applications,  and  others  are  suitable  for  more  general  use.  A  contact  is 
defined  by  identifying  (via  parts,  part  sets,  segment  sets,  and/or  node  sets)  what  locations  are  to 
be  checked  for  potential  penetration  of  a  slave  node  through  a  master  segment.  A  search  for 
penetrations,  using  any  number  of  different  algorithms,  is  made  every  time  step.  Though 
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sometimes  it  is  convenient  and  effective  to  define  a  single  contact  that  will  handle  any  potential 

contact  situation  in  a  model,  it  is  permissible  to  define  any  number  of  contacts  in  a  single  model 
(Hallquist,  2006). 


In  a  typical  FSI  problem,  the  deformations  can  be  very  large  and  predetermination  of  the  contact 
is  essential.  Contact  options  that  are  available  are:  (i)  One-way  treatment  of  contact  (ii)  Two-way 
treatment  of  contact  (iii)  Tied  contact  (translational  and  rotational  DOF,  with  and  without  failure, 
with  and  without  offset)  (iv)  Single  surface. 


There  are  several  contact-related  parameters  in  that  can  be  used  to  modify  or,  in  many  cases, 
improve  contact  behavior.  Some  of  the  contact  parameters  that  can  be  used  pertaining  to  the  test 
cases  are  (i)  Thickness  offsets,  SLTHK  (ii)  Shell  thickness  offset  recommendations  (iii)  Contact 
sliding  friction,  FS  and  FD  (iv)  Penalty  scale  factors,  SFS  and  SFM  (v)  Contact  thickness  SST 
and  MST  (vi)  Viscous  damping  and  most  importantly  (vii)  Bucket-sort  frequency,  BSORT 
which  refers  to  a  very  effective  method  of  contact  searching  method  to  identify  potential  master 
contact  segments  for  any  given  slave  nodes.  The  Bucket-sorting  is  an  expensive  part  of  the 
contact  algorithm  so  the  number  of  bucket  sorts  should  be  kept  to  a  minimum  to  reduce  runtime. 
BSORT  specifies  the  number  of  time  steps  between  bucket  sorts.  If  two  relatively  smooth 
simply-connected  surfaces  are  moving  across  each  other  without  folds,  the  bucket  sorting  has  to 
be  done  at  larger  intervals.  In  certain  contact  scenarios  where  contacting  parts  are  moving 
relative  to  each  other  in  a  rapid  fashion,  such  as  airbag  deployment,  more  frequent  (than  default) 
bucket  sorting  intervals  improves  the  contact  and  impact  behavior  (Hallquist,  2006). 


The  scaled  model  for  the  wave  basin-actual  dimensions  of  the  movable  beach  and  the  membrane 
structure  (Figure  6)  shows  various  components  of  the  numerical  wave  basin  including  a 
wavemaker,  movable  beach  and  MOSES  (membrane  structure  containing  internal  fluid  -  water 
in  this  case).  The  water  domain  and  the  air  domain  form  an  essential  component  of  a  typical  FSI 
model.  The  air  domain  (with  zero  density)  is  needed  to  capture  the  free  surface  of  the  water. 
Figure  7  shows  an  isometric  view  of  the  numerical  wave  basin  and  Figure  8  shows  a  close-up 
view  of  the  movable  beach  and  the  membrane  structure  set-up. 
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Part  3:  Wavemaker  Part  2:  Air  Domain 


r  Beach  (Solid  Wedge)  structure  and  springs 

L  x  attached  to  membrane 

structure 


Figure  6  Components  of  numerical  wave  basin 
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Figure  7  Isometric  view  of  numerical  wave  basin 
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Figure  8  Close-up  view  of  movable  beach  and  membrane  structure 

The  different  parts  that  are  used  in  modeling  the  numerical  wave  basin  are  listed  below.  The 
units  of  measurement  used  to  build  the  geometry  and  for  the  analysis  are  N,  m,  kg  and  sec  (SI). 

1.  Water  Domain: 

Material:  Null  Hydrodynamics 
Section:  Solid  elements  (SOLIDALE) 

Equation  of  State:  Gruniesen  EOS 

2.  Air  Domain: 

Material:  Vacuum  (with  zero  density) 

Section:  Solid  elements  (SOLID  ALE) 

3.  Wavemaker: 

Material:  Rigid  material 
Section:  Shell  elements 

4.  Movable  Beach: 

Material:  Rigid  material 
Section:  Solid  elements 

5.  Membrane  Structure: 

Material:  Rigid  material  and  Flexible  (Material  Elastic) 

Section:  Shell  elements 

6.  Springs: 

Material:  Spring  elastic 
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Section:  Discrete  elements 


LS-DYNA  SPECIFIC  MODELING  PROCEDURE 


The  pre-processing  procedure  involves  nine  key  steps  demonstrating  sequentially  the 
construction  of  the  modeling  of  the  numerical  wave  basin.  The  nine  steps  are: 

STEP  1:  Modeling  Wave  Basin  and  Various  Components 

STEP  2:  Material  Models 

STEP  3:  Element/Section  Properties 

STEP  4:  Assigning  Material  and  Section  properties 

STEP  5:  Boundary  Conditions 

STEP  6:  Contact  Definition 

STEP  7:  Defining  Curves  for  Wavemaker  Motion  and  Gravity 

STEP  8:  Application  of  ALE  Multimaterials 

STEP  9:  Control  Cards  for  ALE  with  Time  Step  allocation 

Note:  A  detailed  description  of  the  sequential  steps  along  with  material  and  section  properties 
will  be  discussed  later  in  the  Appendix  section. 


PRELIMINARY  SIMULATION  RESULTS 


Numerical  simulations  were  conducted  for  15sec  (model  time).  The  models  were  tested  with 
available  material  characterization  for  the  air-water  domain  and  the  membrane  structure.  Results 
shown  here  are  the  sequential  evolutions  of  the  modeling  of  the  wave  basin.  The  results  have 
been  classified  into  four  different  cases.  Each  case  is  designed  to  validate  specific  features  of  the 
numerical  model  with  increasing  level  of  complexity  in  order  to  validate  each  modeling 
component  to  capture  the  real  physics  of  the  experiment. 

Case-I  shows  the  behavior  of  the  membrane  structure  under  the  action  of  waves  generated  in  the 
NWB  with  a  reflective  wave.  Case-II  shows  the  improvements  made  to  Case-I  by  introducing  a 
non-reflective  BC  and  added  friction  between  the  membrane  structure  and  the  beach. 
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Plots  shown  here  are  the  d3plot  output  files  captures  using  a  screen  dump  from  LS-PrePost.  The 

time  taken  to  run  Case-I,  II  and  III  varied  from  12  to  48  hours  on  10  CPUs  on  the  Oregon  State 
Universities  College  of  Engineering  clusters. 


An  isolated  case  of  the  membrane  structure  treated  as  flexible  is  shown  in  Case-IV.  The 
termination  time  for  this  case  was  lOsec  and  the  job  ran  on  a  Naval  computer  network  on  40 
nodes  (CPUs)  for  approximately  4  days. 


Case-I 


The  membrane  structure  in  this  case  is  treated  as  rigid  (K=5000N/m  with  added  weight  to  the 
nodes  where  the  springs  are  located).  The  stroke  of  the  wavemaker  is  0.1524m.  The  wave 
reflects  back  from  the  far  end  of  the  wave  basin.  Friction  between  the  membrane  structure  and 
the  beach  was  neglected.  Figure  45  shows  the  wave  propagation  above  the  membrane  structure 
in  the  wave  basin. 
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Figure  9  Wave  propagation  above  the  rigid  membrane  structure  (Case-I) 


Case-II 

In  this  case,  the  membrane  structure  is  also  treated  as  rigid  (K=5000N/m  with  added  weight  to 
the  nodes  where  the  springs  are  located).  The  stroke  of  the  wavemaker  is  0.1524m.  The  above 
case  had  a  reflecting  boundary  but  the  present  case  has  a  non-reflecting  boundary  with  friction 
considered  between  the  beach  and  membrane  structure.  As  expected,  the  wave  reflects  back  from 
the  far  end  of  the  wave  basin,  which  confirms  the  reflecting  feature  of  Case-I.  Figure  46  shows 
the  wave  propagation  above  the  membrane  structure  for  the  present  case. 
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Figure  10  Wave  propagation  above  the  rigid  membrane  structure  (Case-II) 
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Case-Ill:  Test  Cases  with  Real  Shape  of  the  membrane  structure 


The  membrane  structure  in  this  case  has  the  actual  shape  and  dimensions  used  in  the  experiment 
(treated  as  a  rigid  body).  Importantly  for  the  present  case  the  wavemaker  has  a  stroke  of  0.0762m 
and  a  frequency  of  1Hz.  The  spring  constant  is  same  as  used  for  the  previous  working  cases 
(K=10,000N/m  with  added  weight  to  the  nodes  10kg  each).  Figure  47  shows  the  wave 
propagation  above  the  membrane  structure  for  Case-Ill.  Note  the  wave  propagation  at  the  left 
edge  of  the  wavemaker/water  domain. 
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Figure  11  Wave  propagation  above  the  rigid  membrane  structure  (Case-Ill) 
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Case-IV :  Flexible  test  case 

The  membrane  structure  in  this  case  is  treated  as  fully  elastic  (Material  Elastic)  with  shell 
elements  with  the  rest  of  the  parameters  remaining  same  as  in  previous  cases.  The  von-Mises 
Stress  contours  on  the  membrane  structure  are  shown  at  different  times  of  the  simulation.  Figure 
48  shows  the  stress  contours  plotted  using  the  LS-PrePost  on  the  flexible  membrane  structure. 
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Figure  12  LS-PrePost  stress  contour  plots  for  a  flexible  membrane  structure  (Case-IV) 
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CONCLUDING  REMARKS 


The  major  objective  of  this  study  is  to  assess  the  feasibility  of  using  an  advanced  numerical  code 
to  simulate  the  dynamic  motion  of  MOSES  under  wave  action.  This  is  achieved  by  developing  a 
numerical  model  of  the  wave  basin  with  wave  generation  over  a  rigid  basin  bottom  to  study  the 
dynamic  behavior  of  MOSES  under  laboratory  conditions.  The  numerical  predictions  will  be 
compared  qualitatively  with  experimental  results  obtained  in  an  accompanying  experimental 
study  by  Bloxom  et  al  2010.  The  simulations  thus  far  utilize  an  Arbitrary  Lagrangian  and 
Eulerian  (ALE)  technique  to  capture  the  multi-physics  phenomenon. 


Numerical  predictions  are  first  used  to  complement  the  original  experimental  data  and  then  the 
results  are  used  to  supplement  experimental  tests  to  establish  trends  for  different  orientations  of 
the  MOSES.  The  first  step  in  the  modeling  was  to  treat  MOSES  as  a  fully  rigid  object  and  tests 
its  motion  under  the  action  of  the  waves.  Results  show  that  the  coupled  FSI  numerical  code  can 
successfully  capture  the  motion  of  MOSES  under  the  severe  action  of  waves. 


As  a  first  step  the  object  was  placed  on  top  of  the  movable  beach  and  it  was  allowed  to  move 
freely  (at  this  point  of  model  testing  the  geometry  of  the  object  (MOSES)  was  not  complete 
owing  to  the  challenge  of  modeling  the  membrane  structure  shape  for  the  object).  It  was  noticed 
that  when  the  object  is  left  to  move  freely  on  the  beach,  it  initially  slides  on  the  beach  and  sinks 
eventually.  This  was,  however,  an  important  step  to  comprehend  the  behavior  of  each  component 
of  the  coupled  system  under  the  action  of  the  waves  in  the  NWB.  The  next  step  in  the  procedure 
was  to  model  and  attach  springs  at  the  end  of  the  object  to  control  the  motion  of  the  inflatable 
causeway  under  the  action  of  waves. 


Once  the  desirable  motion  of  the  object  is  achieved,  the  subsequent  step  was  to  change  the 
stiffness  and  the  weight  of  the  object  to  obtain  a  more  realistic  behavior  of  the  motion.  The 
waves  generated  in  the  numerical  wave  basin  were  allowed  to  reflect  back  after  they  reach  the 
end  of  the  NWB  to  observe  the  motion  of  the  object  under  the  reflective  wave. 


The  numerical  test  cases  were  made  even  more  realistic  by  introducing  friction  between  the 
beach  and  the  object  to  capture  the  physics  correctly.  Once  the  spring/weight  were  successfully 
modeled  and  the  motion  established  with  an  added  friction  the  following  simulations  were 
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performed  under  the  realistic  case  of  a  non-reflecting  wave  motion  replicating  a  real  ocean 

scenario  and  the  results  show  accurate  physics  of  the  motion  of  MOSES. 


All  the  test  cases  described  above  show  good  matching  of  the  physics  of  the  wave  motion  and 
the  motion  of  the  membrane  structure  on  the  movable  beach  and  the  qualitative  behavior 
compared  well  with  those  shown  in  the  experimental  test  cases.  Both  the  numerical  test  cases  of 
treating  the  water-fdled  membrane  (MOSES)  as  rigid  and  fully  flexible  show  a  realistic  behavior 
of  the  membrane  structure  motion  under  the  action  of  waves.  However,  the  case  where  the 
membrane  structure  was  treated  as  fully  flexible  was  computationally  very  intensive  and  it  is 
hard  to  achieve  a  number  of  cycles  of  the  wave  motion. 


Given  the  limited  time  frame  available  to  develop  a  numerical  model  and  conduct  simulation  for 
the  inflatable  causeway,  the  present  study  was  carried  out  using  existing  computational  support 
available  at  the  time.  The  computing  needed  to  complete  the  numerical  simulations  ranged  from 
3  hours  (with  8  CPUs  on  the  OSU  cluster)  to  3  days  (20  CPUs  on  a  Naval  distributed  computer 
cluster).  In  the  future,  if  and  when  a  robust  high  performance  computing  platform  is  available  the 
numerical  analysis  of  the  inflatable  causeway  can  be  tested  for  a  fully  flexible  case  of  the 
membrane  structure.  Specifically,  the  study  can  be  further  extended  to  the  cases  of  varying 
orientation  of  the  membrane  structure  and  the  movable  beach. 


31 


Center  for  Innovation  in  Ship  Design 
A  Feasibility  Study  on  Numerical  Modeling  of  Large 
Scale  Naval  Fluid  Structure:  Contact  Impact  Problems 

REFERENCES 


Bloxom,  A.,  Medeling,  A.,  Vince,  C  and  Yim,  C.S.,  2010.  “Modeling  and  testing  of  inflatable 
structures  for  rapidly  deployable  port  infrastructures,”  Technical  report  NSWCCD-CISD- 
2010/010-Ship  System  Integration  and  Design  Department,  Naval  Surface  Warfare  Center, 
Carderock  Division,  West  Bethesda,  MD  20817-5700. 

Belytschko,  T.,  Flanagran,  D.  F.,  and  Kennedy,  J.  M.,  1982.  “Finite  element  method  with  user- 
controlled  meshes  for  fluid-structure  interactions,”  Computer  Methods  in  Applied  Mechanics  and 
Engineering,  33,  pp.  689-723. 

Dickens,  K.,  and  Rosen,  P.,  2007.  “MOSES-Inflatable  Causeway”  Technical  report  NSWCCD- 
CISD-2007 '/005-Ship  System  Integration  and  Design  Department,  Naval  Surface  Warfare  Center, 
Carderock  Division,  West  Bethesda,  MD  20817-5700. 

Hallquist,  J.  O.,  2006.  “LS-DYNA  Theory  Manual,”  Livermore  Software  Technology 
Corporation. 

Tokura,  S.,  and  Tetsuli,  I.,  2005.  “Simulation  of  wave  dissipation  mechanism  on  submerged 
structure  using  fluid-structure  coupling  capability  in  LS-DYNA”,  5th  European  LS-DYNA  users 
conference,  pp.  2c-37. 

Souli,  M.,  Olovsson,  L.,  and  Do,  I.,  2002.  “ALE  and  fluid-structure  interaction  capabilities  in 
LS-DYNA,”  7th  International  Users  Conference,  May  19-21,  Dearborn,  Michigan,  USA. 


32 


Center  f on  Innovation  in  SHip  Design 


Knowledge 


People 


Innouation 


Center  for  Innovation  in  Ship  Design 
A  Feasibility  Study  on  Numerical  Modeling  of  Large 
Scale  Naval  Fluid  Structure:  Contact  Impact  Problems 


APPENDIX 

LS-DYNA  SPECIFIC  MODELING  PROCEDURE 


In  this  appendix,  details  of  the  numerical  simulation  procedure  specific  to  modeling  the  physical 
wave  basin  and  water-filled  membrane  (MOSES)  using  a  multi-physics  multi-scale  finite- 
element  code  LS-DYNA  are  delineated.  Specifically,  a  nine-step  simulation  procedure  pertaining 
to  such  class  of  fluid-structure  interaction  problems  are  shown. 


STEP  1:  Modeling  Wave  Basin  and  Various  Components 


The  wave  basin  consists  of  the  following  domains: 

1)  Water  Domain 

2)  Air  Domain 

3)  Wavemaker 

4)  Movable  Beach 

5)  Membrane  Structure 

6)  Springs 

Creating  Objects  in  LS-DYNA 

Dimensions  of  Water  and  Air  Domain  of  Numerical  Wave  Basin  (NWB): 

The  original  length  of  the  water  domain  is  13.8m  (45.5  ft)  and  the  width  of  the  water  domain  is 
5.6m  (18.5  ft).  The  height  of  the  water  column  is  0.254m  (10  in).  Figure  9  shows  the  elevation  of 
the  water  domain  in  2-D.  The  length  of  the  water  domain  has  been  discretized  into  139  Elements 
(with  an  element  size  of  0.1m).  The  wavemaker  was  located  2m  (82in)  from  the  left  end  of  the 
wave  basin.  The  height  of  the  wavemaker  is  0.7112m  (28in)  and  its  width  is  that  of  the  wave 
basin.  Isometric  view  of  the  numerical  wave  basin  is  shown  in  Figure  10. 
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Figure  13  Elevation  of  Waver-Air  domain  of  actual  wave  basin 
(Green:  Water  domain  and  Red:  Air  domain) 


Figure  14  Isometric  view  of  actual  wave  basin 


Numerical  Model  of  Wave  Basin: 

In  the  physical  experiment,  there  was  material  present  at  the  back  of  the  wave  basin  to  provide 
energy  absorption  to  prevent  wave  reflection.  The  physical  nature  of  the  material  is  not  important 
and  the  absorption  effect  can  be  modeled  numerically  using  an  absorbing  boundary  condition. 
The  run  time  of  the  numerical  model  can  be  significantly  reduced  without  sacrificing  accuracy  in 
capturing  the  physics  of  the  wave  motion  and  the  motion  of  the  membrane  structure  by  placing 
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the  absorbing  boundary  sufficiently  farther  downstream  of  the  structure.  In  this  study,  a  water 

domain  length  of  4.98m  (16.3ft)  is  deemed  sufficient.  Water  and  air  domain  has  been  created 
with  solid  elements.  Elevation  of  the  water  domain  is  shown  in  Figure  11. 


Eight-node  brick  elements  and  four-node  Belytschko-Tsey  shell  elements  have  been  used  for 
discretization  of  the  NWB.  Solid  brick  elements  are  used  for  air/water/movable  beach  domains 
and  shell  elements  for  membrane  structure/wavemaker  domains. 


Part  1:  Water  Domain: 

The  dimensions  of  the  wave  basin  are  as  follows: 

Length  of  the  water  domain  =  4.98m  (16.3ft) 

Width  of  the  water  domain=  5.6m  (18.5  ft) 

Height  of  the  water  domain=  0.254m  (10  in) 

Size  of  the  element=0.1m  (so  the  no.  of  elements  per  length  is  50) 
The  elevation  of  the  water  domain  is  shown  in  Figure  11. 


Figure  15  Elevation  of  water  domain  (2D-view) 
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Part  2:  Air  Domain 

The  dimensions  of  the  air  domain  are  as  follows: 

Length  of  the  air  domain  =4. 98m  (16.3ft) 

Width  of  the  air  domain=5.6m  (18.5  ft) 

Height  of  the  air  domain=0.72m  (30in) 

Size  of  the  element=0.1m  (so  the  no.  of  elements  per  length  is  50) 


Figure  16  Elevation  of  air  domain  (2D-view) 


Note:  In  order  to  visualize/capture  the  free  surface  in  LS-DYNA  an  air  domain  on  top  of  the 
water  domain  is  required.  A  criteria  to  choose  the  height  of  the  air  domain  is  based  on  the 
maximum  wave  height  (usually  the  height  of  the  air  domain  is  chosen  in  such  a  way  that  it  is 
larger  than  the  maximum  wave  height).  Figure  12  shows  the  elevation  of  the  air  domain  in  2-D 
and  Figure  13  shows  the  elevation  of  both  the  air  and  water  domain  in  the  numerical  wave  basin. 
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Visualization  of  Water  and  Air  domain 


Figure  17  Elevation  of  air  and  water  domain  (2D-view) 

Figures  14  and  15  show  three  dimensional  views  of  the  water  and  air  domains.  Figure  15  shows  a 
transparent  view  of  the  air  and  water  domian.  This  feature  is  a  key  in  visualizing  the  wave 
motion  and  the  movement  of  the  membrane  structure  once  the  wavemaker  starts  moving.  It  can 
also  be  used  to  plot  the  velocity  vectors  across  the  domain. 


Figure  18  Isometric  view  of  air  and  water  domains  (3D) 
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Figure  19  Transparent  view  of  air  and  water  domains 
Part  3:  Wavemaker 

A  piston-type  wavemaker  is  used  to  generate  prescribed  regular  waves  or  irregular  waves  from 
specific  spectra.  The  wavemaker  has  a  6”  maximum  stroke  and  the  frequency  of  the  waves  has  a 
2  Hz  maximum.  For  the  current  simulation  a  stroke  of  3”  with  a  maximum  frequency  of  1Hz  is 
used.  The  wavemaker  is  made  of  shell  elements  (Figure  16).  The  Belytschko-Tsay  formulation 
(Hallquist,  2006)  is  used  for  the  shells  and  the  thickness  of  the  shell  is  0.05m.  The  most 
important  part  is  the  orientation  of  the  normals  of  the  wavemaker.  Specifically,  the  normals  have 
to  be  in  the  direction  of  the  wave  generation  (Figure  17).  The  height  of  the  wavemaker  is 
0.71 12m  (28in)  and  the  model  width  is  that  of  the  physical  wave  basin. 
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Figure  20  Isometric  view  of  face  of  wavemaker  and  front  view  of  water  domain 


Figure  21  Vector  normals  of  wavemaker  in  line  with  water  domain 
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Part  4:  Movable  beach 


The  movable  beach  is  made  of  solid  elements  (Figure  18)  and  has  the  following  dimensions: 

Height=0.24m  (9.75”) 

Width  =2.4m  (96”) 

Length  =2. 08m  (82”) 

Slope=l:8.4 


In  this  preliminary  study,  the  beach  is  modeled  as  fixed  to  capture  the  physics  of  the  motion  of 
the  membrane  structure  correctly. 
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Figure  22  Movable  Beach 

Part  5:  Membrane  Structure 

The  membrane  structure,  made  up  of  shell  elements,  is  placed  on  top  of  the  movable  beach. 
Modeling  of  the  membrane  structure  is  one  of  the  most  challenging  parts  of  the  current  study  (in 
addition  to  achieving  the  correct  physics).  The  dimensions  of  the  membrane  structure  are  3.048m 
(10ft)  in  length,  0.2286m  (0.75ft)  in  width,  and  0.2286m  (0.75ft)  in  height,  respectively.  Figure 
19  shows  the  membrane  structure  on  top  of  the  beach. 
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Figure  23  Membrane  structure  and  movable  beach 

After  all  the  five  parts  have  been  modeled  independently  they  are  then  assembled  to  get  the 
components  ready  for  numerical  simulation.  Figure  20  shows  the  assembly  of  air  and  water 
domain,  wavemaker,  solid  wedge  (movable  beach)  and  the  membrane  structure.  Figure  20  shows 
the  side  view  of  the  assembly.  An  isometric  view  of  the  complete  assembly  is  shown  in  Figure 
21. 


Figure  24  Complete  assembly  of  air,  water,  wavemaker,  beach  and  membrane  structure 
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Figure  25  Up-close  of  complete  assembly  (view  from  the  wavemaker  side) 


Part  6:  Constraint  Spring  Modeling 

Two  springs  were  attached  to  the  membrane  structure  to  restrain  one  end.  The  idea  of  modeling 
springs  is  essential  in  the  current  simulation  as  without  them  the  membrane  structure  would  slide 
freely  on  the  moveable  beach  and  move  with  the  wave  which  is  not  desirable  and  practical.  With 
the  springs,  one  end  of  the  membrane  structure  is  restrained.  However,  the  membrane  structure  is 
free  to  rotate  about  the  z-axis  along  the  restrained  end. 

Modeling  Technique: 

The  springs  have  been  modeled  using  discrete  elements  (by  picking  four  duplicate  end  nodes  of 
the  membrane  structure).  These  nodes  become  the  basis  for  the  formation  of  the  spring  elements 
and  they  have  been  assigned  a  mass  value  of  10kg  each.  Once  the  springs  are  modeled,  the  next 
step  is  to  detach  two  nodes  and  extend  them  in  the  z-direction  to  set  the  spring  form.  The 
remaining  two  free  nodes  of  the  spring  are  then  merged  with  the  membrane  structure.  The  spring 
constant  (K)  of  10,000N/m  seem  to  hold  the  membrane  structure  in  place  (a  trial-and-error 
approach  was  used  in  determining  the  spring  constants)!  Figure  22  shows  an  up  close  view  of  the 
springs  and  the  membrane  structure.  It  is  important  to  note  that  the  springs  are  treated  as  fully 
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flexible  to  allow  them  to  move  freely  with  the  motion  of  the  bag.  A  isometric  view  of  the 

movable  beach  with  the  springs  attached  to  the  membrane  structure  is  shown  in  Figure  23. 


As  a  part  of  the  membrane  structure  rests  on  the  bottom  of  the  water  basin  friction  has  to  be 
defined  between  the  membrane  structure  and  the  water  basin  and  it  has  to  be  defined  along  with 
the  springs  attached  to  the  end  of  the  membrane. 


Figure  26  Up-close  of  springs  attached  to  membrane  structure 
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Figure  27  Isometric  view  of  movable  beach  and  springs  attached  to  membrane  structure 

Finally,  a  snapshot  of  the  complete  model  information  with  the  total  number  of  nodes,  elements 
and  materials  is  shown  in  Figure  24. 


Figure  28  Complete  model  information  for  simulation  of  numerical  wave  basin 
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STEP  2:  Material  Models 


Part  1:  Water  Domain 

As  recommended  by  the  LS-DYNA  theory  manual,  to  model  fluid  domain  “Material_Null”  is 
used.  For  solid  elements  the  “Material  Null”  is  the  material  model  that  has  to  be  used  to  invoke 
the  equations  of  state  (EOS)  to  avoid  deviatoric  stress  calculations.  (Hallquist,  2006).  The  Null 
material  model  demands  the  use  of  Gruniesen  equation  of  state  for  a  complete  definition  of  fluid 
contact  modeling.  The  properties  of  water  that  are  needed  for  this  card  are  as  follows: 

Density  =1000  kg/m3 

Pressure  cut-off  —  1.0E+26  N/m2 

Viscosity=1.0E-03  N*s/m2 

As  mentioned  above  the  water  domain  needs  an  equation  of  state  to  capture  the  fluid  behavior 
(which  is  implemented  by  defining  an  EOS  card).  In  physics  and  thermodynamics,  an  EOS  is  a 
relation  between  state  variables.  It  is  a  constitutive  equation  which  provides  a  mathematical 
relationship  between  two  or  more  state  functions  associated  with  the  matter,  such  as  its 
temperature,  pressure,  volume,  or  internal  energy.  Equations  of  state  are  useful  in  describing  the 
properties  of  fluids,  mixtures  of  fluids,  solids,  and  even  the  interior  of  stars. 

Detailed  description  of  the  equation  of  state  for  water  (the  Gruniesen  equation)  is  described  as 
follows. 

The  Gruniesen  equation  of  state  with  cubic  shock  velocity-particle  velocity  defines  pressure  for 
compressible  material  is  governed  by 


p0C2jU 

1  + 

fi-A 

2  J 

a  2 

m-2m\ 

1-fo  -1)/*-' 

c  A2  c  A3 

2  A  + 1  3(a  +  02J 

+  (r0+av)E 


(1) 


where  E  is  the  internal  energy  per  initial  volume,  C  (speed  of  sound)  is  the  intercept  of  the  shock 
velocity-particle  velocity  curve  (Us-Up),  S) ,  S2  and  S}  are  the  coefficients  of  the  slope  of  the 
us-up  curve,  y0  is  the  Gruniesen  gamma,  and  a  is  the  first  order  volume  correction  to  y0 .  The 
constants  C,  S, ,  S2,  S3,  y0  ,  and  a  are  all  input  parameters.  Compression  is  defined  in  terms  of 

the  relative  volume,  V,  as:  p  =  --!■  For  expanded  materials,  the  pressure  is  defined  by: 
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P  =  p{)C~ /.i  +  (y  (l  +aju)E .  p{)  is  the  reference  density.  The  input  used  in  the  Gruniesen  equation  is 
given  in  Table  1  (Hallquist,  2006). 


Table  1:  Gruniesen  Parameters  for  Water 


Equation  of  State 

Card  value 

(Parameter) 

(from  manual) 

C 

100  m/s 

S, 

1.979 

s2 

0.0 

s3 

0.0 

r0 

0.11 

A 

0.0 

E0 

0.0 

K 

0.0 

Part  2:  Air  Domain 

Material_Vacuum  is  used  for  the  air  domain.  The  only  parameter  this  card  demands  is  the 
density  of  air  and  the  value  is  set  to  zero. 


Part  3:  Wavemaker 

The  wavemaker  is  treated  as  rigid.  Hence,  Material_Rigid  is  used  for  this  purpose.  The 
properties  of  steel  are  generally  used  and  Figure  25  shows  the  value  inputs  for  this  card.  It  is 
important  to  mention  that  CMO,  CONI  and  CON2  are  the  values  which  assign  constraints  on  the 
global  motion  of  the  wavemaker.  The  values  of  1,  5  and  7  specify  that  the  motion  of  the 
wavemaker  is  only  in  the  positive  x-direction  interacting  with  the  fluid. 

Note: 

>  Use  “CMO”  (Center  of  Mass  Constraint  Option)  to  set  BC’s  of  both  the  mound  and 
wavemaker 

>  “CON  1”  =  Translational  constraint 

>  “CON2”  =  Rotational  constraint 

>  The  wavemaker  uses  “5”  and  “7”  (allows  it  to  move  in  only  X-translational  direction 
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Figure  29  Material  properties  for  wavemaker 

Part  4:  Movable  Beach 

The  movable  beach  is  made  up  of  solid  elements  and  it  is  treated  as  rigid.  Hence,  just  like  the 
wavemaker  the  beach  uses  the  Material_Rigid  with  the  same  properties  as  of  the  wavemaker. 
The  only  major  difference  however  is  that  the  beach  is  fixed  globally  so  the  values  of  CMO, 
CONI  and  CON2  are  1,7  and  7,  respectively. 

Part  5:  Membrane  Structure 

The  membrane  structure  too  is  treated  as  rigid.  So  Material  Rigid  is  used  with  CMO,  CON  1  and 
CON2  as  1,  0  and  0  respectively  as  the  membrane  structure  is  free  to  move.  The  properties  of 
steel  used  for  the  wavemaker  and  movable  beach  are  retained. 

Part  6:  Springs 

The  springs  are  modeled  first  as  rigid  and  then  as  a  flexible  part.  For  rigid  case  as  mentioned 
above  Material_Rigid  is  used  with  the  properties  of  steel  and  when  it  is  treated  as  a  flexible  part 
Material_Spring_Elastic  is  used.  Figure  26  shows  the  card  values  for  the  elastic  case. 
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Figure  30  Material  properties  for  springs 
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STEP  3:  Element/Section  Properties 


Eight-node  brick  elements  and  four-node  Belytschko-Tsey  shell  elements  have  been  used  for 
discretization  of  the  NWB.  Solid  brick  elements  are  used  for  air/water/movable  beach  domains 
and  shell  elements  for  membrane  structure/wavemaker  domains.  The  Belytschko-Lin-T sey  shell 
element  is  implemented  in  LS-DYNA  as  a  computationally  efficient  alternative  to  the  Hughes- 
Liu  shell  element  (Hallquist,  2006)  used  earlier.  For  a  shell  element  with  five  through-the- 
thickness  integration  points,  the  Belytschko-Lin-Tsey  shell  elements  requires  725  mathematical 
operations  compared  to  4066  operations  for  the  under-integrated  Huges-Liu  shell  elements. 
Because  of  its  computational  efficiency,  the  Belytschko-Lin-Tsey  shell  element  is  usually  the 
shell  element  formulation  of  choice.  The  Belytschko-Lin-T  sey  shell  element  is  based  on  a 
combined  co-rotational  and  velocity-strain  formulation.  The  efficiency  of  the  element  is  obtained 
from  the  mathematical  simplifications  that  result  from  these  two  kinematical  assumptions.  The 
co-rotational  portion  of  formulation  avoids  the  complexities  of  non-linear  mechanics  by 
embedding  a  coordinate  system  in  the  element.  The  choice  of  velocity-strain  rate  of  deformation 
in  the  formulation  facilitates  the  constitutive  evaluation,  since  the  conjugate  stress  is  the  more 
familiar  Cauchy  stress.  Detailed  formulation  of  this  element  is  given  in  LS-DYNA  theory  and 
keyword  manual. 

Part  1:  Water  Domain 

The  water  domain  is  modeled  using  an  ALE  formulation  which  forms  the  basis  for  all  FSI 
problems.  Figure  27  shows  the  card  details  for  Solid  ALE. 


Figure  31  Section  properties  for  water  domain 
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Part  2:  Air  Domain 

The  air  domain  is  also  modeled  using  SolidALE.  The  card  values  as  shown  in  Figure  25  are 
retained  for  the  air  domain.  However,  the  air  domain  does  not  need  an  equation  of  state. 


Part  3:  Wavemaker 

The  wavemaker  is  modeled  using  shell  elements.  Figure  28  shows  the  card  values  for 
Section_Shell.  Belytschko  Tsay  element  formulation  (ELFORM-2)  is  used  for  this  purpose  and 
the  thickness  of  the  shell  is  chosen  as  0.05m. 


Figure  32  Section  properties  for  wavemaker 


Part  4:  Movable  Beach 

The  beach  is  modeled  using  regular  solid  elements.  The  following  are  some  of  the  options  that 
are  available  for  the  solid  element  (Hallquist,  2006).  Figure  29  shows  the  section  properties  for 
the  beach. 

ELFORM:  Element  formulation  options 

EQ.O:  1  point  corotational  for  *MAT_MODIFIED_HONEYCOMB 

EQ.l:  constant  stress  solid  element  (default) 

EQ.2:  fully  integrated  S/R  solid 

EQ.3:  fully  integrated  quadratic  8  node  element  with  nodal  rotations 
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EQ.4:  S/R  quadratic  tetrahedron  element  with  nodal  rotations 
EQ.5:  1  point  ALE 
EQ.6:  1  point  Eulerian 

EQ.7:  1  point  Eulerian  ambient,  EQ.8:  acoustic 

EQ.9:  1  point  corotational  for  *MAT_MODIFIED_HONEYCOMB 

EQ.10:  1  point  tetrahedron  etc. 

EQ.  1  i.e.,  the  Constant  Stress  Solid  Element  is  used  in  the  present  case. 
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Figure  33  Section  properties  for  movable  beach 


Part  5:  Membrane  Structure 

The  membrane  structure  is  modeled  using  shell  elements  and  the  same  formulation  as  used  for 
the  wavemaker  is  used.  The  membrane  weight  in  the  present  study  is  negligible  compared  to  the 
weight  of  water. 


Part  6:  Springs 

Springs  are  modeled  using  discrete  elements.  Figure  30  shows  the  card  details  for  the  discrete 
elements.  The  default  of  zero  are  used  for  all  the  card  values. 
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Figure  34  Section  properties  for  defining  spring  elements 
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STEP  4:  Assigning  Material  and  Section  properties 


Figure  31  shows  the  way  to  assign  the  various  material  and  section  properties  to  various 
components.  This  forms  a  very  important  part  of  the  modeling  as  a  proper  way  of  assigning  these 
cards  can  lead  to  error  free  runs. 


Figure  35  Assigning  material/section/EOS  data  to  all  the  parts  in  LS-PREPOST 
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The  boundary  conditions  of  the  numerical  model  can  be  specified  by  using  *Set_Nodes.  To 
define  the  boundary  conditions  for  the  fluid  domain  in  a  convenient  way,  the  *Set_Node  option 
is  used  to  combine  nodes  that  have  the  same  BC’s.  First  highlight  only  those  parts  to  which  the 
modeler  wants  to  create  a  node  set.  Figure  32  shows  the  boundary  conditions  defined  for  wave 
basin.  The  figure  highlights  the  constraints  applied  to  the  nodes  on  the  wall  of  the  wave  basin 
only. 


1 

2 

3 

4 

5 

6 

7 

8 


Fix_XYZ  (nodenum=8)(sub:l) 
Fix_XY  (nodenum=220)(sub:l 
Fix_XZ  (nodenum=36)(sub:l) 
Fix_YZ  (nodenum=196)(sub:l 
Fix_X  (nodenum=495)(sub:l) 
Fix_Y  (nodenum=5390)(sub:l 
Fix_Z  (nodenum=882)(sub:l) 
Fix_XYZ_Spring  (nodenum=2) 


These  are  the  eight  node  sets 
created  to  define  the  constraints  on 

the  nodes  at  various  locations  in 

the  wave  basin 
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Figure  36  Application  of  boundary  conditions 
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Boundary_SPC_Set  Definitions 


Once  the  node  sets  have  been  created  using  the  node  set  option,  the  next  step  is  to  assign  these 
nodes  to  the  Boundary_SPC_Set  card  to  complete  the  process  of  defining  various  constrains  of 
motion  on  the  nodes  of  the  wave  basin.  Figure  33  shows  the  SPC  set  card.  In  the  figure,  the 
Fix  XYZ  node  constraints  are  shown. 


Figure  37  Boundary  SPC  set  card  definition 


Non-Reflecting  Boundary  Condition  Definition 

The  end  of  the  wave  basin  where  there  is  no  wavemaker,  the  boundary  is  specified  as  non¬ 
reflecting.  In  order  to  define  the  non-reflection  boundary  condition,  the  first  and  foremost  task  is 
to  create  a  segment  set  for  the  non-reflective  side.  Figure  34  shows  the  segment  sets  that  have 
been  created.  Figure  35  shows  the  card  where  the  non-reflection  BC  is  defined. 
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Figure  38  Segment  set  created  to  define  the  non-reflecting  BC 


Figure  39  Non-reflecting  BC  card 
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STEP  6:  Contact  Definition 


Contact  definitions  form  an  integral  part  of  a  typical  FSI  problem.  In  the  present  study,  the 
“Constrained-Lagrange-In-Solids”  type  of  Interface/contact  is  used  to  model  the  impact  event 
between  the  beach/membrane/wavemaker  (Lagrangian)  and  the  water-air  part  sets  (Eulerian). 
Penalty  coupling  without  erosion  is  used  in  the  present  numerical  model  because  the  fluid 
elements  are  solids  and  structural  ones  are  shells. 

Constrained-Lagrange-In-Solids  couples  a  Lagrangian  mesh  (slave)  of  shells  to  the  material 
points  of  an  Eulerian  mesh  (master  surface).  In  this,  the  moving  surface  of  three-dimensional 
beach/membrane/wavemaker  (a  Lagrangian  mesh)  is  treated  as  the  slave  surface,  and  the  target 
water-air  mesh  is  treated  as  the  master  surface.  The  slave  part  or  slave  part  set  is  coupled  to  the 
master  part  or  master  part  set.  The  recommended  input  penalty  factor  is  0.1.  The  quadrature  rule 
for  coupling  slaves  to  master  employs  a  rectangular  grid  of  4  x  4  points.  The  coupling  direction  is 
normal  to  the  structure  (wavemaker)  and  the  stress  at  the  interface  is  compression  only. 

In  addition,  all  the  fluid  elements  are  set  up  to  be  an  ALE  multi-material  group  so  that  each 
element  has  a  specific  material  property  and  the  ALE  formulation  tracks  the  motion  of  the 
elements  in  the  numerical  simulation.  Note  that  the  ALE  control  settings  are  essential.  In  our 
numerical  model,  the  number  of  time  steps  between  mesh-smoothing  and  advection  cycles, 
which  could  speed  up  run  time  but  reduce  stability,  is  determined  to  be  unity  based  on  the 
consideration  of  computational  stability.  The  first  order  advection  method  is  applied  and  it  is 
efficient  and  accurate  for  the  water.  An  option  is  selected  (the  AFAC  value  is  set  at  -1)  to  ensure 
that  no  significant  element  distortion  is  allowed  in  the  numerical  simulation;  only  expansion  and 
contraction  are  permitted  (LS-DYNA  theory  manual). 

Hence,  there  is  contact  between: 

Set  1 :  Wavemaker  and  Beach  (master  set)  and  Air  and  Water  domains  (slave  set) 

Set  2:  Membrane  structure  (master  set)  and  Air  and  Water  domains  (slave  set) 

Figure  36  shows  the  card  details  of  the  contact  between  the  two  sets. 
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Figure  40  Constrained  Lagrangian  in  Solid  card  details 

It  is  important  to  note  that  the  NQUAD  value  is  set  to  4  which  determine  the  number  of  coupling 
points  assigned  to  each  surface  of  the  Lagrangian  elements.  In  addition,  the  normal  direction  of 
the  wavemaker  shell  elements  is  very  important  for  the  coupling  to  take  place  (Figure  37). 


Figure  41  NQUAD/Normal  direction  of  coupling  between  Lagrangian  and  Eulerian  parts 


57 


Center  f on  Innovation  in  SHip  Design 


Knowledge 


People 


Innouation 


Center  for  Innovation  in  Ship  Design 
A  Feasibility  Study  on  Numerical  Modeling  of  Large 
Scale  Naval  Fluid  Structure:  Contact  Impact  Problems 

STEP  7:  Defining  Curves  for  Wavemaker  Motion  and  Gravity 


>  Application  of  Gravity  Globally 

The  gravitation  force  is  applied  quasi-statically  over  a  2-second  duration  prior  to  the  dynamic 
simulation  and  it  remains  constant  throughout  the  simulation.  The  value  of  gravity  acceleration  is 
9.81m/s2.  Figure  38  shows  the  graph  for  the  application  of  gravity.  It  is  however  important  to 
mention  that  in  LS-DYNA  quasi-static  application  of  gravity  is  recommended  to  capture  accurate 
and  stable  pressure  measurements  in  the  water  domain  (this  is  the  main  reason  for  the  variation 
of  gravity  from  t=0  to  t=2sec  in  the  graph  shown). 
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Figure  42  Gravity  loading  on  wave  basin  (Acceleration  time  history) 


>  Wavemaker  Motion 

A  piston  type  wavemaker  is  used  in  the  current  simulation.  The  stroke  of  the  wavemaker  is  set  at 
3”  with  a  frequency  of  1Hz.  The  waves  generated  are  regular  waves  with  simple  irregular  wave 
spectra.  Figure  39  shows  the  graph  for  the  wavemaker  motion  displacement  with  time.  It  is 
important  to  note  that  the  wavemaker  motion  starts  after  gravity  has  set  in  properly  (i.e.,  after  2 
seconds). 
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Figure  43  Wavemaker  displacement  time  history 


Once  the  curves  for  gravity  and  wavemaker  are  defined  the  next  step  is  to  assign  them 
independently.  The  gravity  curve  is  assigned  by  using  the  Load  Body  Y  option  and  the 
wavemaker  motion  is  assigned  by  used  the  Boundary  Prescribed  Motion  Rigid  card.  Figure  40 
shows  the  Load  Body  card.  This  card  is  used  to  set  the  gravity  force.  The  Boundary  Prescribed 
Motion  Rigid  card  is  used  to  prescribe  the  wavemaker  motion  (Figure  41).  Note  that  the 
wavemaker  moves  only  in  the  x-translational  direction.  VAD  is  set  to  2  which  defines  the  motion 
in  the  form  of  displacement. 


Figure  44  Load  Body  card  to  define  gravity  loading 
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■  KEYWORD  INPUT 


Figure  45  Boundary  prescribed  motion  rigid  card  for  wavemaker  motion 


>  Defining  Friction  Between  Beach  and  Membrane  Structure 

Figure  42  shows  the  define  friction  card  between  the  beach  and  the  membrane  structure.  As  a 
part  of  the  membrane  structure  rests  on  the  bottom  of  the  water  basin  friction  has  to  be  defined 
between  the  membrane  structure  and  the  water  basin. 
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KEYWORD  INPUT 


Figure  46  Friction  between  membrane  structure/beach  and  bottom  of  wave  basin 
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In  this  model,  both  the  water  and  air  domains  are  treated  as  ALE  parts  so  they  form  an  ALE 
multimaterial  set.  Figure  43  shows  the  ALE  multimaterial  card  set  definition.  In  the  card  shown, 
“1”  stands  for  the  water  domain  and  “2”  stands  for  the  air  domain. 


Figure  47  Multimaterial  ALE  card 
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STEP  9:  Control  Cards  for  ALE  with  Time  Step  allocation 


Control  ALE  is  an  important  card  that  defines  the  type  of  continuum  treatment  needed  for  the 
ALE  formulation  along  with  the  method  of  mesh  advection.  This  card  has  to  be  defined 
whenever  and  wherever  an  ALE  formulation  is  used  in  the  numerical  simulation.  Figure  44 
shows  the  card  showing  various  card  values  for  the  control  ALE  feature.  The  termination  time  is 
set  to  10sec-25sec  depending  on  the  problem  complexity  and  the  time-step  is  set  to  0.9. 


Figure  48  Control  ALE  card  details 


Details  of  some  of  the  card  parameters  are  tabulated  below. 


DCT:  Default  continuum  treatment: 
EQ.  1 :  Lagrangian  (default) 

EQ.2:  Eulerian 

EQ.3:  Arbitrary  Lagrangian  Eulerian 
EQ.4:  Eulerian  Ambient 


METH:  Advection  method 

EQ.  1 :  donor  cell  +  HIS  (first  order  accurate) 

EQ.2:  Van  Leer  +  half  index  shift  (second  order) 

EQ.-2  Modified  Van  Leer 

EQ.3:  donor  cell  +  HIS,  first  order  accurate, 
conserving  total  energy  over  each  advection  step 
instead  of  conserving  internal  energy 


The  time  interval  between  the  outputs  is  defined  by  using  the  *  Database  card.  The  outputs  are  in 
the  d3plot  format  and  the  “dt”  value  is  set  to  0.05  for  the  current  simulation.  Finally,  check  for 
the  duplicate  nodes  and  perform  a  reference  check  and  see  if  there  are  any  unassigned  card 
values.  Ultimately  renumber  the  node/element/material  and  section  numbers  before  saving  the 
file. 
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